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Abstract. Several attempts have been made to detect solar-like oscillations in the GO IV star rj Boo. We present 
here new observations on this star simultaneously conducted with two spectrographs: Coralie mounted on the 
1.2-m Swiss telescope at the ESO La Silla Observatory (Chile) and Elodie based on the 1.93-m telescope at 
the Observatoire de Haute Provence (France). In total, 1239 spectra were collected over 13 nights. The power 
spectrum of the high precision velocity time series clearly presents several identifiable peaks between 0.4 and 
1.0 mHz showing regularity with a large and small separations of Av = 39.9 /xHz and 81/02 ~ 3.95 /iHz respectively. 
Twenty-two individual frequencies have been identified. Detailed models based on these measurements and non- 
asteroseismic observables were computed using the Geneva evolution code including shellular rotation and atomic 
diffusion. By combining these seismological data with non-asteroseismic observations, we determine the following 
global parameters for 77 Boo: a mass of 1.57 ± 0.07 M0, an age t = 2.67 ± 0.10 Gyr and an initial metallicity 
{Z/X)i — 0.0391 ± 0.0070. We also show that the mass of 77 Boo is very sensitive to the choice of the observed 
metallicity, while the age of rj Boo depends on the input physics used. Indeed, a higher metallicity favours a higher 
mass, while non-rotating models without overshooting predict a smaller age. 

Key words. Stars: individual: rj Boo - Stars: oscillations - techniques: radial velocities 



1. Introduction 

The measurements of the frequencies of p-mode oscilla- 
tions provide an insight into the internal structure of stars 
and are nowadays the most powerful constraint to the the- 
ory of stellar evolution. The five-minute oscillations in the 
Sun have led to a wealth of information about the solar 
interior. These results stimulated various attempts to de- 
tect a similar signal on other solar-like stars by photomet- 
ric or equivalent width measurements, with little success 
due to the extreme weakness of the expected amplitude. 
These past years, the stabilized spectrographs developed 
for extra-solar planet detection achieved accuracies needed 
for solar-like oscillation detection by means of radial veloc- 
ity measurements (Carrier et al. 120031 Bouchy & Carrier 

A primary target for the search for p-mode oscilla- 
tions is the well-studied bright subgiant GO 7] Bootis 
(HR5235). Several attempts have been made to detect 
solar-like oscillations in the GO IV star r] Boo. The first 
result was obtained by Kjeldsen et al. H1995|l with observa- 
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tions conducted with the 2.5-m Nordic Optical Telescope 
(NOT) on La Palma. In contrast to all other detections 
which were based on velocity measurements obtained us- 
ing high-dispersion spectrographs with stable references, 
they monitored changes in the equivalent widths (EW) 
of temperature-sensitive spectral lines. This enabled them 
to determine a large separation of 40.3 fiHz. Meanwhile, a 
search for velocity oscillations in 77 Boo using the AFOE 
spectrograph by Brown et al. (1997) has failed to detect 
a stellar signal. The analysis of NOT data was refined by 
Kjeldsen et al. I|2()()3(l using all existing complementary 
data: new EW measurements obtained with the NOT, 
new radial velocities (RV) measured at Lick Observatory 
and RVs from Brown with the AFOE spectrograph. They 
found a large separation of — 40.4 /iHz and identified 
21 oscillation frequencies. 

In this paper, we report Doppler observations of r/ 
Bootis made with the Coralie and the Elodie spectro- 
graphs in a multi-sites configuration. These new measure- 
ments confirm the detection of p-modes and enable the 
identification of twenty- two individual mode frequencies, 
which are compared with those independently identified 
by Kjeldsen et al. (2p03_). Wc also present new models of 
77 Boo based on our seismological constraints. The obser- 
vations and data reduction are presented in Sect. [21 the 
acoustic spectrum analysis and the mode identification in 
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Fig. 1. Radial- velocity measurements of rj Boo. The dis- 
persion reaches 3.8 ms^^. Coralie and Elodie data are 
represented with full circles and open triangles respec- 
tively. The upper panel represents the photon noise un- 
certainties. 



Table 1. Distribution and dispersion of Doppler measure- 
ments. Indications for Elodie measurements are in brack- 
ets. 



Date 


No spectra 


No hours 


a (ms '') 


2002/04/23 


96 


7.83 


3.69 


2002/04/24 


98 


8.00 


3.87 


2002/04/25 


86 


7.96 


3.17 


2002/04/26 


93 


7.54 


3.60 


2002/04/27 


64 


5.23 


3.34 


2002/04/28 


51 


4.35 


3.61 


2002/04/29 


46 (62) 


11.81 


4.52 (4.81) 


2002/04/30 


93 


7.59 


2.97 


2002/05/01 


99 


8.02 


3.86 


2002/05/02 


99 


7.83 


3.56 


2002/05/03 








2002/05/04 


68 (61) 


12.61 


3.98 (4.21) 


2002/05/05 


71 (61) 


12.26 


3.79 (4.01) 


2002/05/06 


91 


7.77 


4.04 



Sect.O the calibration and modeling of f] Boo in Sect.01 
and the conclusion is given in Sect.[Sl 



2. Observations and data reduction 

J] Boo was observed in May 2002 simultaneously with the 
spectrographs Coralie at La Silla Observatory (Chile) 
and Elodie at the Observatoire de Haute Provence 
(France) in order to improve the window function and to 
make the mode identification easier (see Sect.Ol. 

2.1. Coralie measurements 

r/ Boo was observed over fourteen nights (April 23 - 
May 07 2002) with Coralie, the high-resolution (50 000) 
echelle spectrograph mounted on the 1.2-m Swiss telescope 
at La Silla, known for the p-mode identification in the 
a Cen system (Bouchy & Carrier l^0U2l Carrier & Bourban 
I2UU3|I . During the stellar exposures, the spectrum of a 
thorium lamp carried by a second fiber is simultaneously 
recorded in order to monitor the spectrograph's stability 
and thus to obtain high- precision velocity measurements. 
A description of the spectrograph and the data reduction 
process is presented in Carrier et al. (2001") and Bouchy 
et al. H200ip . Exposure times of 180 s, thus cycles of 295 s 
with a dead-time of 115 s, allowed us to obtain 1055 spec- 
tra, with a typical signal-to-noise ratio {S / N) in the range 
of 115-260 at 550 nm. For each night, radial velocities were 
computed relative to the highest signal-to-noise ratio op- 
tical spectrum obtained in the middle of the night. The 
mean for each night is then subtracted. The radial velocity 
measurements are shown in Fig. ^ and their distribution 
and dispersion are listed in Table The dispersion of 
these measurements reaches 3.6ms^^. 



2.2. Elodie measurements 

Due to bad weather, only 184 spectra were collected over 
3 nights with the high resolution (42 000) spectrograph 
Elodie (Baranne et al. I199f)|l mounted on the 1.93-m 
telescope at Haute-Provence Observatory (France). The 
observations are also achieved in simultaneous-Thorium 
mode. The wavelength coverage of the spectra is 3890- 
6815 A, recorded on 67 orders. The dead-time is about 
120 s as the exposure time varies between 150 and 240 s 
depending on the seeing, the extinction and the airmass 
in order to reach a S/N in the range 150-300 at 550 nm. 
The radial velocity determination method is the same than 
for Coralie data. The radial velocity measurements are 
shown in Fig. ^ and their distribution and dispersion are 
listed in Table ^ The dispersion of these measurements 
reaches 4.3 ms~^. 

3. Power spectrum analysis 

In order to compute the power spectrum of the velocity 
time series, we use the Lomb-Scargle modified algorithm 
(Lomb [T976I Scarrie ri982(l with a weight being assigned 
to each point according its uncertainty estimate. The time 
scale gives a formal resolution of 0.87/iHz. The resulting 
periodogram, shown in Fig. |21 exhibits a series of peaks 
near 0.8 mHz, exactly where the solar-like oscillations for 
this star are expected. Typically for such a power spec- 
trum, the noise has two components: 

— At high frequencies it is flat, indicative of the Poisson 
statistics of photon noise. The mean white noise 
level cTpow calculated between 1.2 and 1.6 mHz is 
0.0179m^s~^, namely 11.9cms~^ in amplitude. With 
1239 measurements, this high frequency noise corre- 
sponds to (Tflv = ^ iV CTpow /4 — 2.35 ms^^. This 
radial velocity uncertainty, larger than for others stars 
with the spectrograph Coralie (see e.g. Eggenberger 
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Fig. 2. Power spectrum of the radial velocity measure- 
ments of 77 Boo. The window function is shown in the 
inset. The noise is represented by the white line. 



et al. l2(J(J4al or Bouchy & Carrier l!!UU2fl . can be mainly 
explained by the high rotational velocity of 77 Boo 
{v sini = 12.8 km s~^, determined from Coralie spec- 
tra). 

— Towards the lowest frequencies, the power should scale 
inversely with frequency squared as expected for in- 
strumental instabilities. However, the computation of 
the radial velocities introduces a high pass filter. 
Indeed, the radial velocities were computed relative 
to one reference for each night and the average radial 
velocities of the night fixed to zero (see Sect.|21). This 
results in an attenuation of the very low frequencies 
which can be seen on Fig. [5] 

The power spectrum presents an excess in the range 0.4- 
1.1 mHz. The combined noise has a value decreasing from 
0.027 to 0.019 m^s-^ (14.5 to 12.2 cms'^) in the above 
mentioned interval (see Fig.|21l. The noise has been deter- 
mined by fitting a function of the type 1/z^^. Note that 
the filtering induced by the radial velocities computation 
does not influence the frequency of the peaks in the range 
0.4-1.1 mHz, but could slightly change their amplitudes. 
The amplitude of the strongest peaks reaches 79 cms~^, 
corresponding to a signal to noise of 6 (in the amplitude 
spectrum). This amplitude is estimated as the height of 
the peaks in the power spectrum with a quadratic sub- 
traction of the mean noise level. The analysis of the 184 
Elodie spectra in addition to the Coralie spectra allows 
us to diminish the first and second daily aliases (11.57 
and 23.14/iHz) by only 9%. Note that a full coverage 
of Elodie data would have allowed us to diminish daily 
aliases by 33 %. 

3.1. Search for a comb-like pattern 

In solar- like stars, p-mode oscillations of low-degree are ex- 
pected to produce a characteristic comb-like structure in 
the power spectrum with mode frequencies i>n.i reasonably 
well approximated by the asymptotic relation (Tassoul 
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Fig. 3. Autocorrelation of the power spectrum under- 
sampled with a resolution of 1.75 /i-Hz and a threshold 
of 0.2 m^s^^. The large splitting is estimated to about 
40 /xHz. 



mm : 



u^^i « Aiy{n + + l)Do . 



(1) 



Here, Dq, which is equal to ^61^02 if the asymptotic rela- 
tion holds exactly, and e are sensitive to the sound speed 
near the core and to the surface layers respectively. The 
quantum numbers n and I correspond to the radial order 
and the angular degree of the modes, and Aiy and 6i>02 
to the large and small separations. To search for period- 
icity in the power spectrum, an autocorrelation is calcu- 
lated and presented in Fig. |2| As the large spacing is not 
strictly constant, we use an undersampled power spectrum 
with a resolution of 1.75 /iHz. We are thus less sensitive to 
small variations of the large spacing versus the frequency. 
Each peak of the autocorrelation corresponds to a struc- 
ture present in the power spectrum. The two strong peaks 
at low frequency close to 11.5 and 23 /iHz correspond to 
the daily aliases. The most probable large spacing is situ- 
ated near 40/iHz for rj Boo. The result of the autocorrela- 
tion confirms the large spacing of 40.4 /^Hz deduced from 
Kjeldsen et al. 



3.2. Mode identification 

The frequencies were extracted using an iterative algo- 
rithm, which identifies the highest peak between 400 and 
1100 /iHz and subtracts it from the time series. First, we 
iterated this process until all peaks with an amplitude 
higher than 3<7 were removed (see Fig.^J. a represents the 
noise in the amplitude spectrum decreasing from 14.5 to 
12.2 cms~^ in the above mentioned interval (see Sect.Ol). 
Peaks with amplitudes below the 3(t threshold were not 
considered since they were too strongly influenced by noise 
and by interactions between noise and daily aliases. This 
threshold, which ensures that the selected peaks have only 
a small chance to be due to noise, gave a total of twenty- 
three frequencies (see Table EJ- Because of the daily alias 
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Fig. 5. Power spectrum of rj Boo with the twenty- two extracted frequencies indicated by dashed hues. The identification 
of each extracted frequency is given in Table 01 



Table 2. Identification of extracted frequencies. Some fre- 
quencies can be either £ — 1 modes or due to noise. 



Frequency 


Mode ID 


S/N 


[mHz] 






512.2 


£ = 1 


3.2 


544.6 - 11.6 = 533.0 


£^0 


3.0 


550.3 


£ = 1 


3.0 


589.9 


£ = 1 


4.6 


622.2 - 11.6 = 610.6 


£ = Q 


3.6 


614.1 + 11.6 = 625.7 


£ = 1 


3.3 


653.8 + 11.6 = 665.4 


£ = 1 


3.3 


669.9 


e = 1 


3.9 


691.3 


£ = 


4.4 


724.5 


£ = 2 


5.2 


728.3 


noise 


3.4 


729.5 


£ = Q 


4.6 


748.5 


£ = 1 


6.5 


777.2 - 11.6 = 765.6 


1 = 2 


5.0 


781.0 - 11.6 = 769.4 


l = Q 


4.1 


775.8 4- 11.6 = 787.4 


I = 1 


3.3 


805.1 


£ = 2 


3.0 


809.2 


£ = Q 


3.7 


834.5 -^11.6 = 846.1 


t = 2 


4.1 


888.7 


1 = 2 


4.0 


891.6 


e = o 


3.2 


947.6 


£ = 1 


3.1 


960.3 -M1.6 = 971.9 


£ = 


3.7 



of 11.57 /iHz, we cannot know a priori whether the fre- 
quency selected by the algorithm is the right one or an 
alias. We thus considered that the frequencies could be 
shifted by ±11.57 //Hz, and made echelle diagrams for dif- 
ferent large spacings near 40 /iHz until every frequency 
could be identified as an ^ = 0, ^ = 1 or ^ = 2 mode. In 
this way, we found an averaged large spacing of 39.9 /iHz. 
It is difficult to identify £ = 1 modes as they appear to be 
mixed modes. Some identified £ = 1 modes could thus be 




iii 
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m 



800 1000 
frequency l/xHz] 

Fig. 4. Top: Original power spectrum of rj Boo. Bottom: 
Cleaned power spectrum: all peaks listed in Table |21 have 
been removed. The dot-dashed, dotted and dashed lines 
indicate an amplitude of 5 a, 4(t and 3ct, respectively. 
Numerous peaks are still present below 3 cr, since no peaks 
have been cleaned below this threshold. These peaks can 
be due to p-mode oscillations and noise or have artificially 
been added by the extraction algorithm due to the finite 
lifetimes of the modes 



rather due to noise, e.g. peaks at 625.7 and 665.4 /iHz. The 
peak at 728.3 /iHz was attributed to the noise, however it 
could be related to the £ = mode 729.5 /tHz "split" owing 
to its lifetime. 

The echelle diagram showing the twenty-two identified 
modes is shown in Fig. |3 The frequencies of the modes 
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Table 3. Oscillation frequencies (in /iHz). The frequency 
resolution of the time series is 0.87 ^Hz. 



£ = £^ 1 1 = 2 



Ai^e 39.9 (0.1) 39.7 (0.2) 40.9 (0.2) 



are shown in Fig. [S] and are given in Table |31 with radial 
order of each oscillation mode deduced from the asymp- 
totic relation (see Eq. 1) assuming that the parameter e is 
near the solar value (e© ^ 1.5). We can see that the oscil- 
lation modes do not strictly follow the asymptotic relation 
due to mixed £ = 1 modes and a large curvature of others 
modes in the echelle diagram. The average small spacing 
has a value of Si'q2 — 3.95 ± 0.9 /iHz. The large spacing 
is separately determined for each value of i and is given 
in the last line of Table 13 The weighted average of these 
three Ai^i yields the value of — 39.9 ± 0.1 /iHz. 

The twenty-two identified modes are compared to pre- 
viously identified ones by Kjeldsen et al. ipDt)3|i) in Fig. |H1 
Both identifications are rather in good agreement at low 
frequency but present major discrepancies at high fre- 
quency. Although the present data have a higher S/N, 
additional measurements are needed to resolve these am- 
biguities. 

3.3. Oscillation amplitudes 

Concerning the amplitudes of the modes, theoretical com- 
putations predict oscillation amplitudes between 1 and 
1.5 ms^^ for a 1.6 Mq star like rj Boo, with mode lifetimes 
of the order of a few days (Houdek et al. I1999|l . The am- 
plitudes of the highest modes, in the range 55-80 cm s~^, 
are then lower than expected. The observations indicate 
that oscillation amplitudes are typically 2.5-3.5 times so- 
lar. This disagreement can be partly explained by the life- 
times of the modes. Indeed, the oscillation modes have 
finite lifetimes, because they are continuously damped. 
Thus, if the star is observed during a time longer than 
the lifetimes of the modes, the signal is weakened due to 
the damping of the modes and to their re-excitation with 
a random phase. 



1000 









_ 


• 


_ 


- 












- 


■ • u 


_ 
- 




, n U 


_ 
- 




, • n 








_ 




• 


_ 












- 








- 


0* 

























□ •) 








▲ 






A 




• 









, , , , , , 



10 20 30 

u modulo Av [fiUz] 

Fig. 6. Echelle diagram of identified modes (in black) with 
a large separation of 39.9 fiHz. The modes i=2 (■), £—0 
(•), and £=1 (a) are represented with a size proportional 
to their amphtude. Open symbols correspond to modes 
determined by Kjeldsen et al. I)2()()3|l . Both identifications 
are in good agreement at low frequency. However, at high 
frequency their £ = 2 mode frequencies(n) seem to be 
identified in this paper as £ = modes (•), and their 
£ = modes (o) are not present in our data. 

4. Comparison with models 

In order to compare the asteroseismic observations with 
theoretical predictions, stellar models were computed us- 
ing the Geneva evolution code including shellular rotation 
and atomic diffusion (see Eggenberger et al. 2005 for more 
details). We used the OPAL opacities, the NACRE nu- 
clear reaction rates ( Angulo et al. 11999(1 and the standard 
mixing-length formalism for convection. 

4.1. Non-asteroseismic observational constraints 

Following Di Mauro et al. ij^UnSIl (hereafter DM03), the 
luminosity of 77 Boo was deduced from the Hipparcos par- 
allax: L/Lq = 9.02 ± 0.22. Concerning the effective tem- 
perature, we added recent results to the references used 
by DM03 to determine a new average value (see Table 0J). 
As a result, the effective temperature T^g = 6030 ± 90 K 
was adopted. This value is in perfect agreement with the 
effective temperature of 6028 ± 45 K used by DM03, with 
a larger error which seems to us more realistic in view of 
the different values found in the literature. The box in the 
HR diagram for ry Boo delimited by the observed effective 
temperature and luminosity is shown in Fig. [7| 

The metallicity of Boo adopted by DM03 is [Fe/H] = 
0.305 ± 0.051. Compared to different observed metallici- 
ties, this value seems to be quite large. We thus decided 



n = 


11 




512.2 




n = 


12 


533.0 


550.3 




n = 


13 




589.9 




n = 


14 


610.6 


625.7 




n = 


15 




665.4 / 669.9 




n = 


16 


691.3 




724.5 


n = 


17 


729.5 


748.5 


765.6 


n = 


18 


769.4 


787.4 


805.1 


n = 


19 


809.2 




846.1 


n = 


20 






888.7 


n = 


21 


891.6 






n = 


22 




947.6 




n = 


23 


971.9 
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Table 4. Metallicity and temperature determination for 
77 Boo (since 1990). The errors on the selected parameters 
are chosen to encompass all acceptable values (last line) . 



Teff [° K] 


\Fe/H] 

L / J 


References 


6000 


0.16 


Mc William 


6219 


0.30/0.37 


Balachandran 1990] 


6068 


0.19 


Edvardsson et al. .19931 


5943 


0.20 


Gratton et al. 119961 


_ 


0.23 


MisheninajlMn 




0.28 


Fuhrma (Cayrel de StrobelHSHJ 


6003 


0.25 


Genarro et al. 120011 




0.16 


Buzzoni et al. I^IMl 


6000 


0.25 


Feltzing & Gonzalez 


6120 


0.16 


GrayPnUTl 


5964 


0.19 


Mallik et al. 


5957 


0.15 


Nordstrom et al. 120041 


5942 


0.18 


AUende Prieto et al. I^Uni 


6030 ± 90 


0.23 ±0.07 





to adopt a lower average value of [Fe/H] — 0.23 ± 0.07, 
which is determined from the recent measurements listed 
in Table H 

Finally, we used the observed surface velocity of 77 
Boo to constrain the rotational velocity of our models. 
From Coralie spectra, we determined a rotational veloc- 
ity vsmi = 12.8 km s~^. Since the value of the angle i is 
unknown, we assumed that it is close to 90°. Thus our 
models of rj Boo have to reproduce a surface velocity of 
about 13kms~^. 

4.2. Computational method 

Basically, the computation of a stellar model for a 
given star consists in finding the set of stellar modeling 
parameters which best reproduces all observational data 
available for this star. The characteristics of a stellar 
model including the effects of rotation depend on six 
modeling parameters: the mass M of the star, its age 
(t hereafter), the mixing-length parameter a = I /Hp 
for convection, the initial surface velocity VJ and two 
parameters describing the initial chemical composition of 
the star. For these two parameters, we chose the initial 
hydrogen abundance Xi and the initial ratio between the 
mass fraction of heavy elements and hydrogen {Z/X)i. 
Assuming that this ratio is proportional to the abundance 
ratio [Fe/H], we can directly relate (Z/X) to [Fe/H] by 
using the solar value {Z/X)q = 0.0230 given by Grevesse 
& Sauval (1998). Moreover, we fixed the mixing-length 
parameter to its solar calibrated value (a© = 1.75) and 
we assumed the initial hydrogen abundance Xi to be 
Xi — 0.7. As a result, any characteristic ^ of a given 
stellar model has the following formal dependences with 
respect to modeling parameters: A — A{M,t, Vi, {Z/X)i). 



The determination of the set of modeling parameters 
(M, t, Vi, {Z/X)i) leading to the best agreement with the 
observational constraints is made in two steps. First, we 
constructed a grid of models with position in the HR di- 
agram in agreement with the observational values of the 
luminosity and effective temperature (see Fig. [TJ. Note 
that the initial ratio between the mass fraction of heavy 
elements and hydrogen {Z/ X)i is directly constrained by 
the observed surface metallicity [Fe/H], while the initial 
velocity Vi is directly constrained by the observed rota- 
tional velocity. 

For each stellar model of this grid, \ow-l p-mode fre- 
quencies were then calculated using the Aarhus adiabatic 
pulsations package written by J. Christensen-Dalsgaard 
H1997|l . Following our observations, modes £ < 2 with fre- 
quencies between 0.4 and 1.1 mHz were computed and the 
mean large (Az^) and small spacings (151^02) were deter- 
mined. The mean large spacing was computed by consider- 
ing only the radial modes. Once the asteroseismic charac- 
teristics of all models of the grid were determined, we per- 
formed a minimization as in Eggenberger et al. I|2004b|l . 
Thus, two functionals are defined: Xtot ^^'^ xLtcro- The 
X^Qt functional is defined as follows: 

^ ^ / ^thco ^obs \ 

Xtot = X! y ' crC°'^''' j ' 

where the vectors C contains the following observables for 
one star: 

C = (i/LQ,Teff, [Fe/H], Aiy,,5i.o2) • 

The vector C*''°° contains the theoretical values of these 
observables for the model to be tested, while the values 
of C"''^ are those listed above. The vector ctC contains 
the errors on these observations. Note that the observed 
rotational velocity is not included in this minimization, be- 
cause of its large uncertainty resulting from the unknown 
inclination angle i. The Xastero functional is defined as fol- 
lows: 

Xastero — I ^ ) ' 

z— 1 ^ 

(3) 

where a = 0.87 fiHz is the error on the observed frequen- 
cies estimated as the frequency resolution, iV = 22 is the 
number of observed frequencies, and (D^) is the mean 
value of the differences between the theoretical and ob- 
served frequencies : 

1 ^ 

1=1 

The determination of the best set of parameters was based 
on the minimization of the functional defined in equa- 
tion 121) which includes three non-asteroseismic and two 
asteroseismic observational constraints. Once the model 
with the smallest Xtot determined, we refined the grid 
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in the vicinity of this preliminary solution in order to find 
the best solution which minimizes at the same time Xtot 
and xLtcro- 



systematic difference {D,^) between theoretical and ob- 
served frequencies has been taken into account. Following 
Christensen-Dalsgaard et al. H1995fl . the theoretical oscil- 
lation amplitudes are estimated by 




3.76 



Log T.„ 

Fig. 7. Evolutionary tracks in the HR diagram for two 
models of r] Boo. The dot and the square indicate the lo- 
cation of the Ml and M2 model respectively. The error box 
in continuous line indicates the observational constraints 
of L and T^s used in our analysis, while the dotted box 
corresponds to the constraints used by Di Mauro et al. 



4.3. Results 

Using the observational constraints listed in Sect. 14.11 
with the observed frequencies listed in Table 13 we per- 
formed the minimization described above. We found 
the solution M = 1.57 ± O.O7M0, t = 2.67 ± O.lOGyr, 
Vi ^ 90kms-i and (Z/X)-, = 0.0391 ± 0.0070. The posi- 
tion of this model in the HR diagram (denoted model Ml 
in the following) is indicated by a dot in Fig.[7| The charac- 
teristics of this model are given in Table |5l The confidence 
limits of each modeling parameter given in Tableware es- 
timated as the maximum/minimum values which fit the 
observational constraints when the other calibration pa- 
rameters are fixed to their medium value. Note that the ra- 
dius deduced for 77 Boo R = 2.72 Rq is in good agreement 
with the interferometric radius of 2.766±O.O39i?0 and 
2.682±O.O43i?0 determined respectively with the Mark 
III optical interferometer (Mozurkewich et al. I2003|l and 
with the VLTI (Kervella et al. i2Q04j . 

Theoretical p-mode frequencies of the Ml model are 
compared to the observed frequencies by plotting the 
echelle diagram (see Fig. |H| . Note that in this figure the 



where A is the surface amplitude and e the normalized 
energy of the mode. Aq{i') and eo{v) are obtained by in- 
terpolating to frequency v in the results of radial modes. 

The model shows that £ ~ 1 modes deviate from the 
asymptotic relation for p-modes. This is a consequence of 
the avoided crossings (Christensen-Dalsgaard et al. 119951 
Guenther & Demarque,1996: DM03; Guenther l^ITOl here- 
after GU04). This results in frequencies which are shifted 
relative to the frequencies expected for pure p-modes. 
This is particularly true for the £ = 1 modes at low fre- 
quency which strongly deviate from the asymptotic rela- 
tion. The observed frequencies for £ = 1 modes seem to 
deviate from the asymptotic relation, which is in accor- 
dance with theoretical predictions. However, Fig.|S|shows 
that the model results are not able to precisely reproduce 
the individual frequencies of these modes at low frequency. 

FigureElalso shows that the agreement between obser- 
vations and theoretical predictions for modes with £ = 
and £ = 2 is good, except for the two £ = 2 modes with 
the smallest and the largest frequency = 724.5 and 
888.7 /iHz). Note that the £ — 2 modes are also infiucnced 
by the avoided crossings. However, the effects of coupling 
become much weaker for these modes than for modes with 
£ = 1, since p-modes with £ — 2 penetrate less deep in the 
stellar interior. 

The variations of the large and small spacing with the 
frequency are given in Fig. |51 Large spacings for £ — 1 
modes are not plotted in Fig. |51 since these modes devi- 
ate too strongly from the asymptotic behaviour of pure 
p-modes. Table El and Fig. |51 show that the mean large 
spacing of the Ml model is in perfect agreement with the 
observed value. The observed variation of the large spac- 
ing with the frequency is also correctly reproduced by the 
model, except for the large value of the £ ^ 2 point close 
to 900 /iHz. TableEland Fig.^also show that the observed 
small spacings are compatible with theoretical predictions. 
The observed mean small spacing is however slightly larger 
than the theoretical one; this is mainly due to the large 
value of the observed small spacing at 724.5 /zHz. 

We conclude that the observed frequencies listed 
in Table |31 are compatible with theoretical predictions. 
Although the general agreement is satisfactory, we also 
note some discrepancies between observed and predicted 
frequencies, especially for the £ — 1 modes at low fre- 
quency. 



8 



F. Carrier et al.: New seismological results on the GO IV 77 Bootis 



Table 5. Models for r] Boo. The Ml and M2 models in- 
clude rotation and atomic diffusion, while the M3 model is 
a standard model computed with an overshooting parame- 
ter = 0.2. The upper part of the table gives the obser- 
vational constraints used for the calibration. The middle 
part of the table presents the modeling parameters with 
their confidence limits, while the bottom part presents the 
global parameters of the star. 
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0.0570 
39.9 
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4 A. Discussion of the results and comparison 
previous studies 



with 



Fig. 8. Echelle diagram for the Ml model with a large 
spacing Av = 39.9 /iHz. Open symbols refer to theoret- 
ical frequencies, while the filled circles correspond to the 
observed frequencies listed in Table O Open circles are 
used for modes with £ = 0, triangles for £ — 1, squares 
for £ = 2 and pentagons ior £ = 3. The size of the open 
symbols is proportional to the relative surface amplitude 
of the mode (see text). Modes with too small surface am- 
plitudes (e.g. g modes) are not shown on this diagram. 



Detailed studies of r/ Boo based on the asteroseismic ob- 
servations of Kjeldsen et al. H2(J03|I have already been per- 
formed by DM03 and GU04. Compared to these studies, 
we notice that our Ml model has a smaller mass. Indeed, 
DM03 found that the mass of 77 Boo is limited to the range 
M = (1.64 - 1.75) Mq, while GU04 proposed two differ- 
ent solutions: a model with a mass of 1.706 Mq which 
has exhausted its hydrogen core, and another model with 
a mass of 1.884 which is still on the main-sequence, 
but is approaching hydrogen exhaustion. These two au- 
thors used the same non~asteroseismic constraints (see 
Sect. I4.1|l . However, contrary to the analysis by DM03 
and contrary to the present work, GU04 used a calibration 
method (the QDG method) which is not limited to mod- 
els with position in the HR diagram in agreement with 
the observational values of the luminosity and effective 
temperature. This explains why GU04 found another so- 
lution with a mass of 1.884 Mq, while DM03 determined 
a mass between 1.64 and 1.75 Mq. Recently, Di Mauro et 
al. PKgl (hereafter DM04) showed that mam-sequence 
models provide a match to the observed location of 77 Boo 
in the HR diagram when overshooting from the convective 
core (ofov > 0.1) is included in the computation. 

The fact that our Ml model is less massive than the 
solution of about 1.7 Mq found by DM03 and GU04 can ei- 
ther be due to the different observational constraints used 



or to the different input physics of the evolution codes. 
Indeed, our models include shellular rotation and atomic 
diffusion contrary to the models calculated by DM03 and 
GU04. For stars more massive than about 1.4 Mq, it is 
necessary to introduce another transport mechanism, like 
the rotationally induced mixing, in order to counteract 
the effect of atomic diffusion in the external layers. When 
only atomic diffusion is included in a star with a thin con- 
vective envelope, helium and heavy elements are drained 
out of the envelope, resulting in too low surface abun- 
dances which are incompatible with observation. This is 
illustrated on Fig. ^| which shows the helium profile in 
the external layers at different age during the evolution 
on the main-sequence for a model including only atomic 
diffusion and for the Ml model which includes shellular 
rotation and atomic diffusion. Fig. ^| shows that rota- 
tionally induced mixing prevents the helium from being 
drained out of the convective envelope. Indeed, the de- 
crease of the surface helium abundance during the main- 
sequence evolution is found to be very small for models 
including rotation and atomic diffusion. 
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Fig. 9. Large and small spacings versus frequency for the 
Ml model. Dots indicate the observed values of the large 
{I — 0) and small spacings, while squares correspond to 
the large spacing determined with £ = 2 modes. 
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Fig. 10. Helium abundance profile in the external layers of 
the star at different age during its evolution on the main- 
sequence. The model on top only includes atomic diffusion, 
while the model on bottom includes atomic diffusion and 
shellular rotation. Apart from the inclusion of rotation, 
the two models have been computed with the same initial 
parameters corresponding to the Ml model. 



4.4.1. Effect of the metallicity 

To investigate the effects of non-asteroseismic observa- 
tional constraints on the solution, we decided to redo the 
whole calibration using the non-asteroseismic constraints 
adopted by DM03 and GU04. The metallicity was in- 
creased Z — 0.04 and a temperature of Toff — 6028 ± 
45 K was adopted. Note that we still used our asteroseis- 
mic constraints for this calibration. In this way, we found 
the solution M = 1.69 ± 0.05 Mq, t = 2.65 ± 0.10 Gyr, 
Vi ^ SOkms^i and (Z/X); = 0.0580 ± 0.0070. The posi- 
tion of this model in the HR diagram (noted model M2 
in the following) is denoted by a square in Fig. |7| The 
characteristics of the M2 model are given in Table We 
conclude that the difference in the adopted value for the 
metallicity explains the different mass determined. Indeed, 
Fig.ilshows that the fact that we used T^s = 6030 ± 90 K 
instead of Toff = 6028 ± 45 K has no significant influence 
on the solution since the Ml model is also included in the 
smaller observational box determined by DM03. 

The higher metallicity used by DM03 and GU04 results 
of course from the larger value of the observed [Fe/H]: they 
adopted [Fe/H] = 0.305 ± 0.051, while we fixed [Fe/H] to 
0.23 ± 0.07 for the Ml calibration. However, we notice 
that it also results from the way one relates the observed 
[Fe/H] to the mass fractions Z and X used in the models. 
Indeed, we directly related {Z/X) to [Fe/H] by using the 
solar value {Z/X)q = 0.0230 given by Grevesse & Sauval 
(1998), while DM03 related Z, and not (Z/X), to [Fe/H]. 
As a result, for a same value of [Fe/H] — 0.305, we deter- 
mined Z — 0.032 while DM03 obtained a higher value of 
Z = 0.040 (for X = 0.7). 

We conclude that the derived mass of 77 Boo is very 
sensitive to the choice of the observed metallicity. When 
a metallicity of [Fe/H] = 0.23 ± 0.07 is adopted, a mass 
of 1.57 ± 0.07 M0 is found. When the higher metallicity 
determined by DM03 is used, we obtain a mass of 1.69 ± 
0.05 M0, in perfect agreement with the results of DM03 
and GU04. 



4.4.2. Effect of the input physics 

Contrary to the masses, the ages of the Ml and M2 mod- 
els are very similar (2.67 and 2.65 Gyr respectively) and 
are therefore not very sensitive to a change of metallicity. 
However, this age is larger than the age of 2.393 Gyr ob- 
tained by GU04 for its solution with a mass of 1.706 M©. 
DM03 pointed out that the age of the models depends 
on the inclusion of overshooting: the age is about 2.3- 
2.4 Gyr without overshooting, and between 2.4-2.7 Gyr in 
presence of overshooting. The age of 77 Boo seems there- 
fore to be sensitive to the input physics used. To inves- 
tigate these effects on the solution, we decided to calcu- 
late models without rotation and atomic diffusion using 
the same observational constraints as the M2 model. In 
this way, we find the solution AI = 1.70 ± 0.05 A/q and 
t = 2.39 ± 0.10 Gyr, in perfect accordance with the re- 
sults of GU04 and DM03. Rotating models predict a larger 
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Fig. 11. Ratio of the mass of the convective core to the 
total mass of the star (^cc) as a function of the central 
hydrogen abundance (Xc), for a standard model without 
overshooting, for a standard model with aov = 0.2 and for 
a model with rotation and atomic diffusion. Apart from 
the inclusion of rotation and overshooting, the three mod- 
els have been computed with the same initial parameters 
corresponding to the M2 model. 



age for ?/ Boo than non-rotating ones. This illustrates the 
fact that, for small initial velocities, rotational effects are 
found to mimic the effects due to an overshoot of the con- 
vective core into the radiative zone. Indeed, the lifetimes 
of rotating models are enhanced with respect to those of 
standard models, because the mixing feeds the core with 
fresh hydrogen fuel. As a result, the exhaustion of hydro- 
gen in the central region is delayed and the time spent on 
the main-sequence increases. This can be seen on Fig. 1111 
which shows the ratio of the mass of the convective core to 
the total mass of the star (^cc) as a function of the central 
hydrogen abundance (Xc). We see that the rotating model 
exhibits a larger convective core for a given Qcc, i.e. for a 
given evolutionary stage on the main-sequence, than the 
standard model without overshooting. In the same way, 
the non-rotating model with aov — 0.2 also exhibits a 
larger convective core on the main-sequence than stan- 
dard models without overshooting. This explains why the 
inclusion of rotation or overshooting increases the lifetimes 
of the model on the main-sequence and hence the deduced 
age for rj Boo. 

Finally, we investigated the solution of a model which 
is still on the main-sequence. As found by DM04, these 
models do not provide a match to the observed Tcff and L 
of rj Boo unless overshooting is included. Thus, our anal- 
ysis using the input physics described above leads to only 
one solution which is in accordance with asteroseismic and 
non-asteroseismic observables: the Ml model which is in 
the post-main-sequence phase of evolution. Using the ob- 
servational constraints listed in Sect. 14. ll we tried to deter- 
mine a model of rj Boo which is still on the main-sequence 
by computing non-rotating stellar models including over- 
shooting. In this way, we found that a model computed 
with an overshooting parameter aov — 0.2 and a mass of 
1.7 Mq enables to match the location of ?/ Boo in the HR 
diagram (see Fig.ll2|). As discussed above, the mass of this 
model (denoted model M3 in the following) is lower than 
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Fig. 12. Evolutionary tracks in the HR diagram for two 
models of 1.7 Mq computed with a different amount of 
overshooting. The dot indicates the location of the M3 
model. The error box in continuous line corresponds to the 
observational constraints of L and TcS used in our anal- 
ysis, while the dotted box corresponds to the constraints 
used by DM03. 



the mass of 1.884 determined by GU04 and the values 
of M = (1.75 - 1.90) Mq found by DM04, because of the 
smaller metallicity used in our analysis. The characteris- 
tics of the M3 model are given in Table [3 
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Fig. 13. Echelle diagram for the M3 model with a large 
spacing Aiy = 39.9 /zHz. Open symbols refer to theoretical 
frequencies, while the filled circles correspond to the ob- 
served frequencies listed in Table |31 Open circles are used 
for modes with £ = 0, triangles for ^ = 1, squares for £ = 2 
and pentagons for £ = 3. 

As already pointed out by GU04 and DM04, the fun- 
damental seismic difference between post-main-sequence 
(Ml) and main-sequence (M3) models concerns the 
avoided crossings. Models in the post-main-sequence 
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Fig. 14. Small separations versus frequency for the M3 
model. Dots indicate the observed values of the small sep- 
arations. 



phase show £ = 1 modes that deviate from the asymptotic 
relation, while models still on the main-sequence show no 
occurrence of avoided crossing. Indeed, for these models 
only modes with a radial order lower than the observed 
ones are mixed (see Fig. |H1 and I13f) . Since observation 
show £ — 1 modes that deviate from the asymptotic rela- 
tion, we conclude that models in the post-main-sequence 
phase of evolution are in better agreement with the as- 
teroseismic measurements than the main-sequence mod- 
els. Moreover, the small separation of the M3 model is 
smaller than the one of the Ml model and is therefore in 
slightly less good agreement with the observed frequencies 
(see Fig. I14II . Note that DM04 also found that post-main- 
sequence models are characterized by larger small separa- 
tions than main-sequence models. Although the Ml model 
constitutes the solution which best reproduced all obser- 
vational constraints, the actual precision on the observed 
frequencies does not enable us to definitively reject the 
solution on the main-sequence. 

5. Conclusions 

Our observations of Boo yield a clear detection of p- 
mode oscillations. Several identifiable modes appear in 
the power spectrum between 0.4 and 1.1 /iHz with average 
large and small spacings of 39.9 and 3.95 /iHz respectively 
and a maximal amplitude of 79cms^^. The global results 
are in rather good agreement with the recent work by 
Kjeldsen et al 1)2003(1 who found large and small spacings 
of 40.4 and 3.00 ijMz. However, the comparison of the indi- 
vidual frequencies leads to some discrepancies. Although 
the present data have a higher S/N, additional Doppler 
measurements using spectrographs like Harps (ESO) and 
forthcoming MOST data with a clean window function 
might help resolving these ambiguities. 

We identified 22 mode frequencies which have been 
compared to theoretical models. The combination of non- 
asteroseismic observations now available for 77 Boo with 
the observed p-mode frequencies listed in Table|31leads to 
the following solution: a model in the post-main-sequence 
phase of evolution, with a mass of 1.57 ± 0.07 M©, an age 
t = 2.67 ± 0.10 Gyr and an initial metalhcity {Z/X)i = 



0.0391 ± 0.0070. We also show that the mass of 77 Boo 
is very sensitive to the choice of the observed metalhc- 
ity and that its age depends on the inclusion of rotation 
and atomic diffusion. Indeed, non-rotating models with- 
out overshooting predict a smaller age of 2.39 ± 0.10 Gyr. 
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